Wigner crystal vs. Priedel oscillations in the ID Hubbard model 
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We analyze the fermion density of the Hubbard model in a finite one- dimensional system. A 
Wigner crystal phase should not be expected, since there cannot be any broken symmetry in one 
dimension and moreover the underlying Luttinger liquid theory is scale invariant. Nonetheless, we 
find a relatively sharp crossover even for moderate short range interactions into a definite Wigner 
crystal region. The results are relevant to ultra cold fermionic gases in optical lattices and finite 
quantum wires. 
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Having been predicted in the early days of quantum 
mechanics, the Wigner crystal[H is one of the simplest 
but most dramatic many-body effects: Due to the long 
range repulsive forces, electrons spontaneously form a 
self-organized lattice at low enough densities and tem- 
peratures much different from a free electron gas. Ex- 
perimental verification has been difficult, but very re- 
cently experimental signatures of a Wigner crystal were 
reported in carbon nanotubesQ. Using ultra cold gases 
in optical lattices it is now also possible to produce well- 
controlled correlated fermion systems in restricted di- 
mensions, albeit with short range interactions 0, Q]. 

The Wigner crystal in one dimension has been dis- 
cussed so far mostly in the context of long-range inter- 
actions. In that case it has been predicted by Schulz 
Q that the density-density correlations corresponding to 
an equally spaced inter-particle distance become domi- 
nant. Numerically a crossover to Wigner density waves 
has been observed at strong long-range interactions with 
very few particles [3] • However, for short-range interact- 
ing fermion systems in one dimension possible Wigner 
crystal signatures have not been addressed theoretically 
so far. 

The prototypical model for short range interacting 
fermions is the repulsive (U > 0) one-dimensional (ID) 
Hubbard model 
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Even though many exact results have been derived for 
this model using the Bethe AnsatzQ, the local densi- 
ties in finite chains with open boundaries cannot yet be 
calculated by exact methods. 

Since the on-site scattering is independent of mo- 
mentum, the system becomes effectively scale invari- 
ant at low energies, so that any crossover towards a 
different physical region would be unexpected. More- 
over, a Wigner crystal region should be unstable, since 
2fcF-density Priedel oscillations 0] are always the slow- 
est decaying correlations due to a limited Luttinger 
parameter [f| 0.5 < K c < 1. Nonetheless, sub-dominant 



oscillations at other wave-vectors also exist as has been 
explicitly shown e.g . for the Hubbard model in a finite 
magnetic field (jl. [To}. 

We now study the density distribution in finite Hub- 
bard chains with hard-wall boundaries by a combination 
of bosonization and numerical density matrix renormal- 
ization group (DMRG) calculations. Despite the fact 
that the interactions are short ranged and of moderate 
strength, we find that a distinct Wigner crystal state 
is always stable at low filling. The results show that 
the scale invariance is explicitly broken. The observed 
Wigner crystal state illustrates that in ID even short 
range interactions have an increasing effect with growing 
inter-particle spacing. The observed crossover is related 
to the so-called spin-incoherent Luttinger liquid 11 1 . 

Density oscillations can be calculated exactly for the 
non-interacting model U = in a finite chain with " open" 
boundary conditions ip a .o — ipa,L+i — 0. The summation 
over standing waves sin(nr^-j-x) results in the well-known 
2kp Friedel oscillations [8f 
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where the Fermi wave- vector kp = ^a^jxj is centered 
between the highest occupied and lowest unoccupied level 
in a system with N fermions with spin. Here x is mea- 
sured in units of the lattice spacing. 

Possible Wigner oscillations can only appear at finite 
repulsive interactions U > when fermions with opposite 
spin avoid each other. In the extreme case considered 
by Wigner all fermions would be arranged at maximum 
distances L/N apart. In the ID case, the system remains 
critical but modulations with wave-vector 2irN/L ~ 4fcp 
will become important as we shall see later. 

The Friedel oscillations are also modified by the 
interactions 12]. This can be seen by using the Lut- 
tinger liquid formalism as the low energy effective theory 
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FIG. 1: Local density for U — At and L = 200 for different fill- 
ing showing the Friedel and Wigner crystal oscillations. The 
solid lines correspond to the theoretical prediction in Eq. ([9]). 



FIG. 2: Local density for TV = 20 and L = 200 showing the 
crossover from 2&f to 4&f oscillations with increasing inter- 
action strength U . 



[l3| after linearization around kp of the fermion fields 
ip ax ss e lkFX tpn ja + e~ lkFX tfji J ^. In particular, using 
bosonization one can identify the Friedel oscillations as 
an expectation value of the operator 



O 



LR 



(e l2k ^[jx)^ R Jx)+h.c 



oc sin(2fcpa; + \/2irK c (p c ) cos(v / 27T( y 9 s ). (3) 

Using the standard mode expansion of the s pin field <p s 
and the charge field cf> c , for finite systems 1J, |l5| the ex- 
pectation value is determined to be 

sin (2kpx) 
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up to logarithmic corrections. The interactions change 
the decay rate of the Friedel oscillations compared to 
Eq. ([2]), which appear to be enhanced for repulsive inter- 
actions K c < 1. However, as we will see, the yet unde- 
termined amplitude also strongly depends on interactions 
and filling fraction. 

The derivation of the Wigner oscillations from 
bosonization is more subtle. They arise from interac- 
tions because of the Umklapp term in the Hamiltonian 
density 



cx cos (^Akpx + y/8irK c ip c ^j , 



(5) 



where g% oc U. In first order perturbation theory this 
operator induces a density expectation value (n)u — 
(n) — (n)o relative to the unperturbed case 

(n(,)),oc/^,V O < O l ^yy^i°) (6) 
Jo „ E a - E 



where \a) are all excited states. By using the known 



mode expansion 15] it is possible to calculate the expec- 
tation values for all bosonic excitations la) with the result 
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where g(x 

f% - x). 
cosAkpx/Akpx - -r 
be approximated as[l 
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Using f°° sin(4fc F y)j/~ 2i ^dy 
■ 2K - J- 0{x~ 2 ^"~ l ) the integral 
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The decay rate for the Akp oscillations is faster than 
for the Friedel oscillations in Eq. (TJ| since K c > 0.5 
for the Hubbard model. The linear dependence on 
g^/kp oc UL/N is only accurate to lowest order in per- 
turbation theory, but the typical oscillatory behavior and 
powerlaw will describe the behavior for any U and filling 
N/L. Alternatively, the ad-hoc inclusion of Ou directly 
in the operator expression for the density is also a valid 
approach The explicit derivation from perturbation 
theory above now provides additional information by in- 
dicating an increase of the amplitude with gz/kp, i.e. with 
larger interactions and smaller filling. However, the ex- 
act amplitude cannot be derived from bosonization, so 
that numerical calculations have to be used. 

We have implemented a DMRG algorithm [l(| for the 
model in Eq. (JXJ) in order to calculate the local density 
in finite systems with a given fermion number N . Typ- 
ical densities at various fillings are shown in Fig. [T] for 



3 




FIG. 3: Crossover of the amplitudes in Eq. (J9J (error < 10%) 
for L — 200 from the DMRG data as a function of filling 
and interaction, d) Crossover points (Ai = A2) in the U- 
A/L-plane showing the scaling with UL/N « 30 for L = 200. 
The solid line separates the Wigner and Friedel regions, e) 
Approximate crossover points UL/N as a function of length. 

U = 4t and L = 200, which clearly exhibit the predicted 
oscillations. Figure [5] shows how the local density at a 
given filling of N/L = 0.1 emerges from the slower Friedel 
oscillations to a Wigner crystal pattern as U increases. 

An accurate data analysis is now possible in terms of 
our analytic predictions from Eqs. 10} and (JSJ) 

sin(2fc F :r) „ sin (^-7^1) x ■ 
n(x)=n -A 1 V > Kc+1 TT^7-(9) 

sm (iTT a; )] 2 [ sin {i^i x )_ 

For arbitrary interactions U > and filling N/L the 
Luttinger parameter 0.5 < K c < 1 can be calculated 
exactly @, @, HH, so that there are only two adjustable 
fitting parameters for the amplitudes in the middle of the 
chain Ai and Ai . For the non- interacting case in Eq. ((2|) 
we have A\ = -j^ti and A2 — 0, which we can use as a 
test of the numerical accuracy. The uniform density hq 
is fixed by the requirement that J n(x)dx — N, so that it 
is not an independent fitting parameter, (e.g. no = XTT 
for f7 = 0). 

Figure [1] shows the quality of typical fits to the DMRG 
data. The oscillations in the middle of the chain are very 
well represented by the analytical expression ([9|), while 
there are small deviations near the edges. Deviations 
from Luttinger liquid theory near boundaries have also 




FIG. 4: Schematic occupied states for a) the nearly free case 
and b) for strong interactions ("Wigner crystal state"). 

been observed before in the context of the local density 
of states [lH . In order to determine the asymptotic am- 
plitudes A\ and A2 as accurately as possible we have 
therefore excluded the first few sites near the ends in the 
fits. The fits are sensitive enough to even confirm the ex- 
act values of the wave- vectors 2fcp and 4fcp — lTTj since 
small deviations of order -j^py already would make the 
quality of the fits considerably worse. 

The results for the Friedel amplitude A\ and the 
Wigner amplitude Ai are shown in Fig. [3] for L = 200. 
The amplitudes show a clear crossover from Friedel os- 
cillations to Wigner crystal waves at low filling. Inter- 
estingly, the Friedel oscillations are suppressed exactly 
when the Wigner crystal waves are strong and vice versa. 
Therefore, it is possible to identify two distinct "regions" 
of Wigner crystal and Friedel behavior. There is no true 
phase transition, but we observe distinctly different phys- 
ical behavior in the two regions. 

From the Luttinger liquid theory it is not a priori obvi- 
ous why Friedel and Wigner oscillations cannot be strong 
simultaneously, but this competition can be understood 
in the momentum space picture. The Friedel oscillations 
in Eq. ([2]) are caused by the sum of standing waves in 
the two spin channels up to k-p as indicated in Fig. U-a. 
The Wigner oscillations can be explained by an analo- 
gous picture, if it is assumed that fermions of opposite 
spin are not allowed to occupy the same orbital as shown 
in Fig.|4]-b. The density of the "Wigner crystal state" is 
again given by a sum as in Eq. ^ where the index now 
runs over only one channel up to twice the wave-vector 
2kp. Therefore, in Eq. © we have exactly K c = |, 
A\ = 0, and A2 = > wn i cn is the maximally possi- 

ble amplitude for Wigner oscillations. It is not possible to 
simultaneously have a Friedel and Wigner crystal state, 
leading to a natural competition between 2/cp and 4/cp 
oscillations in agreement with Fig. [3] 

In our case, the Wigner crystal state in Fig. |4]-b still 
has an equal number of spin up and down fermions. This 
is in contrast to charge density waves in ultracold gases 
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and Wigner crystals in higher dimensions, which typi- 
cally tend to be spin polarized. In fact, it appears that 
the Wigner crystal state in Fig. 0]-b can be realized for 
any magnetization, so that the correlations in the spin 
sector must be irrelevant. A small amplitude of spin- 
boson correlations also explains that the expectation of 
the Friedel operator in Eq. ([4]) is strongly suppressed 
in that regime. Hence, the Wigner crystal state in ID 
signals the onset of spin-incoherent Luttinger liquid be- 
havior with a nearly degenerate spin sector and a truly 
different physical behavior (llT|. 

The results in Fig. [3] now indicate when the crossover 
A\ « Ai towards the spin-incoherent Wigner crystal 
state occurs. We observe scaling behavior in the limit 
of small filling in terms of the variable UL/N as shown 
in Fig. [3]-d for L — 200, where the crossover occurs along 
a line with UL/Nt ss 30. This means that a larger aver- 
age particle distance L/N is equivalent to an increase of 
U. At first sight it appears rather counter-intuitive that 
the on-site interaction U should show a stronger effect 
as the average inter-particle distance L/N is increased. 
This behavior is special to one dimension since the total 
kinetic energy scales with (N/L) 3 at low filling, which be- 
comes always smaller than the total interaction energy, 
which scales with (N/L) 2 as N/L — > 0. 

From the Bethe ansatz we find that the Luttinger pa- 
rameter for the Hubbard model also follows a scaling re- 
lation at low filling 



K r w 0.5 



' Nt N 
ML, 



41n2 +0{ 



L 2 ' U 2 - 



(10) 



However, the length dependence of the crossover points 
UL/N in Fig. [5]-e violates the scaling behavior. For a 
given density N / L and interaction strength U all param- 
eters in the Luttinger liquid theory (v c , v s , K c ) are fixed. 
Nonetheless, it is still possible to observe a crossover 
from Wigner to Friedel oscillations as a function of length 
(moving horizontally in Fig. [3}e) . It is quite surprising 
that a critical model can cross over to different physical 
behavior as a function of length only, when all relevant 
parameters are fixed. This remarkable violation of scale 
invariance is not due to higher order operators, but comes 
from the competition between vastly different velocities 
in the spin and the charge sector. Note, that the spin- 
incoherent Luttinger liquid can be defined even at T = 
as the length- dependent region when the ultraviolet spin 
cutoff v s /a becomes smaller than the infrared charge cut- 
off v c /L. In our analysis the curve in Fig. [3]-e only sig- 
nals the first onset of Wigner crystal behavior, so that 
we cannot predict the full analytic behavior from simple 
arguments. In general, it appears that the competition 
of spin and charge energy scales is a common behavior 
in Luttinger liquid physics. Therefore, the crossover be- 
tween the different density oscillations discussed above 
is a generic signature of spin-charge separation in one 
dimension. 



From the experimental side, Luttinger liquid behav- 
ior has so far only been seen in very special cases, 
such as carbon nanotubes 3, 2(| or cleaved edge over- 
growth wires [2l|. There is some hope now that Lut- 
tinger liquid physics can also be realized with ultra-cold 
fcrmionic atoms in nearly ideal geometries formed by op- 
tical traps 22j , which would have the advantage that the 
density distribution discussed here could in principle be 
detected directly using high resolution cameras, electron 
beams [13], or noise interference 24 1. Fermionic gases can 
already be cooled down to less than 1/10 of the Fermi en- 
ergy. The finite temperature will lead to a faster decay 
of the oscillations from the edges that can be accounted 
for in the theory [3]. In fact it would be interesting to 
study the system in a regime where all spin excitations 
are smaller than the temperature. This would be a per- 
fect realization of the spin incoherent Luttinger liquid, 
leading to a complete vanishing of the Friedel oscillations 
while the Wigner oscillations remain. Hard edges can 
also be implemented by focused laser beams or trapped 
impurity atoms. Quantitative details of these issues will 
have to be discussed in future works. 

In summary, we have systematically analyzed the lo- 
cal density distribution as a function of filling, interaction 
strength U, and system size in finite Hubbard chains. A 
combination of bosonization and DMRG calculations al- 
lows a detailed description of the density oscillations in 
terms of the quantitative formula (|9|) . For small interac- 
tions and large fillings 2/cp Friedel oscillations Ai dom- 
inate, while the Wigner crystal amplitude A2 remains 
small. However, for smaller filling or increasing inter- 
actions the overall amplitude A\ of the Friedel oscilla- 
tions is strongly reduced while A2 grows. This signals 
the crossover to a different physical region, which is de- 
scribed by a Wigner crystal state with no double occu- 
pancy of spin up and down fermions. The density oscilla- 
tions we have described here are the most accessible fea- 
ture to study the crossover towards the spin-incoherent 
Luttinger liquid in detail, e.g. using ultra cold fermionic 
gases in ID optical traps. 
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